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I.  Introduction 

This  is  the  final  report  of  work  completed  under  contract 
No.  N00019-76-C-0674.  The  objective  of  this  contract  was  to  de- 
velop the  technique  of  Molecular  Stuffing  to  produce  low-loss 
parabolic  index  optical  fibers.  This  work  was  a continuation  of 
work  done  under  Contracts  Nos.  N00019-76-C-0083  and  N00019-75-C-0074 , 
under  which  much  of  the  theoretical  and  preliminary  experimental 
work  was  done.  The  concluding  work  under  these  contracts  consisted 
of  the  following: 

(i)  Put  into  operation  experimental  profiling  procedures  to 
test  the  theory  developed  under  Contract  N000-9-76-C-0083 . 

(ii)  Institute  a program  for  the  measurement  and  characteri- 
zation of  measured  index  profiles  in  preforms  and  fibers  in  order 
to  determine  the  modifications  needed  in  Part  (i)  above. 

(iii)  Measure  physical  properties  of  cesium  nitrate  solutions 
for  a precise  control  of  unstuffing  bath  concentration  and  its 
variation  in  time. 

(iv)  Measurement  of  attenuation,  dispersion,  bending  and 
scattering  losses  in  our  fibers. 

II  Brief  Review  of  the  Molecular  Stuffing  Process 

Molecular  Stuffing  is  based  on  the  Vycor  process  developed 
more  than  forty  years  ago  at  Corning  but  it  incorporates  a number 
of  important  modifications  that  make  the  process  suitable  for 
optical  fiber  applications.  The  process  utilizes  thermodynamic 
and  chemical  properties  to  produce  a high  purity  glass  preform 
with  a controlled  index  profile.  Preforms  have  been  made  by 
Molecular  Stuffing  that  yielded  fibers  with  optical  attenuation 
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less  than  10  dB/km  and  fibers  having  a graded  index  profile. 

The  basic  steps  in  the  process  are  illustrated  by  the  flow 
diagram  shown  in  Figure  1. 

A.  Glass  Melting 

An  alkali-boros ilicate  glass  is  melted  in  a conventional  re- 
sistance heated  furnace  in  a platinum  crucible.  The  batch  ma- 
terials are  reagent  or  high  purity  commercial  grade  materials 
which  have  transition  metal  impurities  of  less  than  1 ppm  and 
costing  on  the  order  of  $1. 00/lb  in  small  quantities.  The  melting 
and  homogenization  Of  the  glass  can  be  easily  adapted  to  a 24  hour 
schedule.  After  homogenization  is  complete,  the  glass  is  drawn 
into  rods  of  convenient  diameter  and  is  ready  for  further  process- 
ing. 

B.  Phase  Separation 

The  glass  rods  are  placed  in  annealing  furnaces  for  a speci- 
fied time  at  a temperature  below  the  immiscibility  temperature  of 
the  glass,  the  exact  temperature  and  time  being  determined  by 
the  composition  being  used.  For  the  compositions  in  the  range  we 
have  examined,  a temperature  near  500°C  and  a time  of  a few  hours 
is  appropriate.  The  purpose  of  this  heat  treatment  is  to  obtain 
a completely  interconnected  microstructure  having  a characteristic 
size  of  a few  hundred  angstroms  for  the  phases  which  consist  of 
an  alkali  borate  phase  and  one  which* is  more  than  90%  silica,  the 
remainder  being  B2O3.  A fortunate  result  of  the  phase  separation 
is  that  the  ionic  impurities  which  include  the  transition  metals 
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tend  to  segregate  into  the  alkali  borate  phase.  This  facilitates 
their  removal  and  the  purification  accomplished  in  the  next  step 
of  the  process. 

C.  Leaching 

The  alkali-borate  phase  is  susceptible  to  attack  by  weak  acid 
solutions  such  as  HC1  while  the  silica-rich  phase  is  not.  If  the 
proper  interconnected  microstructure  has  been  obtained  in  the  pre- 
vious step,  leaching  for  a time  on  the  order  . of  one  day  will  yield 
a high  silica  porous  structure  free  of  transition  metal  impurities. 
This  porous  structure  is  then  further  cleansed  of  leaching  pro- 
ducts and  acid  by  washing  in  filtered,  deionized  water. 

D.  Dopant  Purification 

Cesium  nitrate  is  available  from  a number  of  suppliers  with 
transition  metal  impurities  ranging  from  tens  of  ppm  to  a few  ppb. 

A number  of  standard  techniques  are  available  for  in-house  purifi- 
cation of  cesium  nitrate  including:  recrystallization,  liquid- 
liquid  extraction,  mercury  hydrolysis,  filtration  and  ion  exchange. 

As  reported  in  the  Quarterly  Report  under  Contract  N00019-76-C-0083 
dated  July,  1976,  both  recrystallization  and/or  liquid-liquid 
extraction  performed  on  Kawecki-Berylco  cesium  nitrate  provides 
sufficient  purification  to  obtain  doped  fibers  with  losses  less 

than  20  dB/km.  Combinations  of  the  above  techniques  are  capable 

0 

of  producing  cesium  nitrate  as  pure  as  that  obtainable  from  Merck 
which  is  in  the  ppb  range  for  both  iron  and  copper.  This  purity 
level  should  be  sufficient  to  yield  fiber  attenuations  below  5 dB/km. 
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E.  Profiling 

The  porous  silica  structure  is  an  ideal  host  for  the  deposi- 
tion of  index  and/or  physical  property  modifying  dopants.  The 
majority  of  our  work  has  been  with  cesium  nitrate  as  a dopant. 

The  dopant  is  dissolved  in  an  aqueous  solution  and  is  allowed  to 
diffuse  into  the  porous  glass  rod.  After  a few  hours  the  dopant 
solution  is  uniformly  distributed  in  the  rod.  A step  index  pro- 
file can  be  obtained  by  first  precipitating  the  cesium  nitrate  by 
means  of  a temperature  drop  and/or  solvent  exchange,  and  then  dis- 
solving the  cesium  nitrate  from  a cladding  layer  in  the  rod.  A 
graded  index  cam  be  obtained  by  replacing  the  dopant  solution 
with  solutions  having  lower  concentrations  of  dopant  and  allowing 
diffusion  to  produce  the  desired  index  profile.  After  the  desired 
profile  is  obtained,  the  dopant  is  precipitated  as  is  done  in  the 
case  of  a step  index. 

P.  Drying  and  Consolidation 

The  stuffed  profiled  rods  are  dried  and  consolidated  in  com- 
bustion tubes  in  which  the  atmosphere  can  be  controlled.  The  rods 
are  slowly  heated  first  to  drive  off  water  vapor,  and  then,  as  the 
temperature  is  raised,  the  cesium  nitrate  decomposes,  giving  off 
nitrous  oxide  and  is  oxidized  to  cesium  oxide..  Final  drying 
and  the  oxidation  of  remaining  iron  impurities  is  accomplished 
by  passing  dry  oxygen  over  the  rods  at  about  l/5th  of  an  atmos- 
phere  pressure.  The  temperature  of  the  rods  is  gradually  raised 
until  the  rod  is  fully  collapsed  at  about  850 °C.  The  rod  preform 
is  then  ready  for  fiber  drawing. 
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Ill  Diffusion  Theory  for  Dopant  Solutions  in  Porous  Preforms 

Under  Contract  N00019-76-C-0083  a theory  was  developed 

for  the  diffusion  of  dopant  solutions  in  a porous  rod.  It  was 

shown  that  by  an  appropriate  choice  of  boundary  conditions  one 

could  produce  a parabolic  dopant  profile  in  porous  glass  rods 

resulting  in  a parabolic  index  profile  in  the  final  glass.  The 

details  of  this  theory  are  given  in  the  final  report  of  the  above 

contract  and  are  summarized  below. 

The  diffusion  equation  to  be  solved  is 

3C  _ r 32C  1 3C 

St  ■ D[sr*-  + r SF1 

where  C = C(r,t)  is  the  dopant  concentration  as  a function  of 
radial  position  and  time  and  D is  the  constant  diffusion  coefficient. 

For  the  cylindrical  geometry  of  a rod  preform  and  boundary 
conditions  consisting  of  a step  change  in  dopant  concentration 
at  the  rod  surface  followed  by  a ramp  decrease  in  concentration 
to  zero  the  solution  to  the  equation  can  be  expressed  as  an 
infinite  series: 

.2  • 


« 8 Jo(Bnx) 

C (r,tf ) = AC  (1-x2 ) + Z 3 T'  » [6n-AC] 
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x = — where  b is  the  radius  of  the  rod; 
b 

Jq  = (gn)  =0,  Jq,  JX  are  Bessel  functions  of  zero  and  first 

order; 
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similar  solutions  are  obtained  for  single  or  multiple  step  changes 


in  concentration  with  time 
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Depending  on  the  choice  of  boundary  conditions  it  has  been 
demonstrated  by  computer  calculation  that  one  can  obtain  a com- 
position profile  whose  standard  deviation  from  a parabola  is 
less  than  10“8  for  the  case  of  a step  change  followed  by  a ramp. 
For  a single  step  the  standard  deviation  can  be  as  small  as  10 
and  for  a double  step  as  small  as  10”3.  mhis  precision  also  ap 
plies  for  profiles  which  represent  small  deviations  from  para- 
bolic which  may  be  required  to  minim izedispersion  in  optical 
waveguides. 


TV.  Profiling  Experiments 

During  this  contract  period,  there  was  a major  effort  to 
produce  the  profiles  predicted  for  the  stuffing  process  by  our 
diffusion  theory.  Our  preliminary  results  showed  that  although 
graded  index  profiles  could  be  achieved,  the  results  were  not  as 
predicted.  It  soon  became  clear  that  variables  not  accounted  for 
in  the  theory  must  be  playing  an  important  role  in  obtaining  the 
final  profile. 

A major  difficulty  in  determining  the  parameters  neglected 
in  the  theory  is  the  fact  that  the  only  step  in  the  process  at 
which  the  dopant  profile  may  be  determined  is  after  the  final  col- 
lapsing. It  is  thus  not  possible  easily  to  determine  at  which 
stage  of  the  process  there  is  a departure  from  our  predictions. . 

The  only  path  open  to  us  was  to  complete  a large  series  of  experi- 
ments systematically  varying  parameters  and  looking  for  an  effect 
on  the  final  product  as  indicated  by  the  index  profile. 

For  reference,  the  details  of  the  basic  stuffing  procedure 
are  outlined  below,  beginning  with  the  leached  porous  rod. 

1.  The  porous  rod  is  soaked  for  three  hours  in  an  aqueous 
solution  of  approximately  66%  cesium  nitrate  by  weight  at  a tem- 
perature near  100 °C. 

2.  Unstuffing  is  accomplished  by  replacing  the  stuffing 
solution  with  98°C  flowing  water  for  7 to  9 minutes  depending  on 
the  diameter  of  the  rod  and  as  required  by  the  diffusion  equations. 

3.  After  unstuffing,  the  flowing  hot  water  is  replaced  by 
flowing  cold  water  (T = 0°C)  for  20  to  30  seconds. 
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4.  Precipitation  is  achieved  by  flushing  the  cold  water  out 
with  methanol  at  3 or  4°C,  followed  by  a soak  for  6 to  15  hours. 

5.  The  rod  is  dried  and  collapsed  as  indicated  earlier. 

The  index  profile  obtained  during  one  of  the  first  tests  of 

this  procedure  is  shown  in  Figure  2.  Although  the  profile  is 
graded,  there  are  several  serious  departures  from  the  predicted 
profile.  The  profile  is  skewed  and  not  centered  on  the  rod  axis. 
The  top  of  the  profile  is  flat  and  the  index  at  the  edge  of  the 
rod  is  much  higher  than  it  should  be.  Tl>is  type  of  profile  was 
typical  of  the  early  results,  although  the  above  mentioned  flaws 
could  be  present  in  varying  intensity.  Some  of  the  resultant  pro- 
files ended  up  more  like  a step  index  than  a graded  index,  as 
seen  in  Figure  3,  for  example. 

The  early  results  seemed  to  indicate  that  some  variable  or 
set  of  variables  present  in  differing  degrees  in  the  rods  seriously 
affected  the  diffusion  process.  It  was  thought  that  something  was 
causing  the  diffusion  coefficient  to  be  smaller  near  the  surface 
of  the  rods  than  it  was  in  the  center.  This  would  explain  the 
high  index  near  the  edge  as  well  as  the  poor  shape  of  the  profile. 
Possible  sources  of  such  a problem  were  identified  and  a series 
of  experiments  performed  to  eliminate  them.  The  experiments  at- 
tempted are  first  identified  and  then  explained  in  detail  below. 

The  rods  were  centerless  ground  after  heat  treatment  to 
remove  any  effect  of  surface  composition  difference  and  to  provide 
a perfectly  cylindrical  geometry. 

2.  The  porous  rods  were  given  a 5%  HF  wash  for  2 to  5 minutes 
to  dissolve  anything  that  might  be  clogging  the  pores  near  the 


surface. 


3.  The  porous  rods  were  washed  with  a .5  N HaOH  solution  for 
15  to  90  minutes  for  the  reason  stated  in  No.  2 above. 

4.  It  was  found  that  silica  gel  formed  in  the  leaching  pro- 
cess could  be  removed  by  washing  the  porous  rods  in  95°C  water  for 
10  days.  It  was  thought  that  the  silica  gel  might  clog  the  pores. 

5.  Several  sources  of  cesium  were  used  for  stuffing.  It  was 
thought  that  some  impurity  in  the  dopant  which  was  not  measured  at 
the  time  ’ (such  as  silica)  could  clog  the  pores. 

6.  The  homogeneity  of  the  base  glass  was  improved  to  elimi- 
nate the  possibility  of  a variation  in  the  micro structure  due  to 
inhomogeneity . 

A.  Centerless  Grinding 

In  the  rod  forming  process,  the  rod  is  drawn  from  the  melt 
with  the  surface  of  the  rod  being  drawn  from  the  top  layer  of 
glass  in  the  melting  pot.  This  surface  layer  can  be  deficient  in 
alkali  and  borate  due  to  volatilization.  This  could  effect  the 
microstructure  near  the  surface  resulting  in  a diffusion  barrier 
there.  A large  number  of  rods  were  centerless  ground  to  remove  the 
surface  layer  and  to  remove  any  effects  that  an  oval  rod  might 
produce.  The  grinding  was  done  after  the  phase  separating  heat 
treatment. 

Aside  from  the  elimination  of  the  obvious  effects  due  to  oval 
rods,  the  grinding  of  the  preforms  seemed  to  have  had  no  effect. 

0 

Figures  3 and  4 show  two  ground  rods  prepared  in  a similar  manner, 
yet  yielding  very  different  results.  The  rod  in  Figure  3 shows  a 
very  flat  top  and  the  index  at  the  edge  is  greater  than  1.47  (it 
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should  be  1.458).  The  rod  in  Figure  4 shows  more  curvature  at 
the  top  and  has  a lower  index  at  the  sides.  These  experiments 
were  repeated  several  times  with  similar  inconclusive  results. 

B.  HF  Wash 

It  was  thought  that  during  the  leaching  or  washing  processes 
something  might  cause  the  surface  pores  to  shrink  or  become  clogged 
with  the  result  that  diffusion  at  the  surface  would  be  restricted. 

A short  (2  to  5 minute)  wash  with  a 5%  Hf  solution  would  allow  the 
pores  near  the  surface  to  be  enlarged  without  affecting  the  central 
pores.  A number  of -rods  were  treated  in  this  way  with  varying  re- 
results. 

One  of  the  best  profiles  that  we  have  obtained  was  made  from 
a rod  having  the  HF  wash.  The  profile  is  shown  in  Figure  5.  A 
least  squares  fit  to  the  central  75%  of  the  rod  shows  a relative 
standard  deviation  from  a parabola  of  about  4%.  A second  rod 
shown  in  Figure  6 was  given  an  identical  treatment.  This  profile 
is  almost  a step  profile.  Although  the  HF  wash  may  help  in  some 
cases,  it  is  not  the  controlling  factor  since  we  were  not  able  to 
reproduce  the  results  shown  in  Figure  5 with  any  consistancy. 

C.  NaOH  Wash 

For  the  same  reasons  as  given  above  for  the  HF  wash,  a 
.5  M NaOH  wash  was  tried.  This  procedure  never  gave  good  results. 

A typical  profile  using  this  procedure  is  shown  in  Figure  7. 

D.  Hot  Water  Wash 

Jit  is  known  that  silica  gel  is  formed  in  the  leaching  pro- 
cess. This  gel  is  deposited  in  the  pores  and  reduces  their  ef- 
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fective  size  thus  affecting  the  diffusion  of  dopant  in  the  rods. 

The  silica  gel  may  be  removed  by  washing  the  rods  for  approximately 
10  days  in  water  at  95°C.  The  removal  of  the  silica  gel  can  be 
observed  by  monitoring  the  weight  of  the  porous  rod.  The  dissolu- 
tion of  the  gel  proceeds  at  a much  greater  rate  than  the  attack 
and  dissolution  of  the  silica  skeleton,  so  it  is  clear  when  remo- 
val is  achieved. 

As  in  the  other  experiments, no  positive  correlation  between 
profile  quality  and  silica  gel  removal  was  seen.  Figure  8 shows 
a typical  profile  achieved.  The  top  is  too  flat  and  the  edges  are 
too  high. 


E.  Dopant  Purity 

The  purity  of  the  dopant,  cesium  nitrate,  in  these  experiment 
varied  considerably.  It  was  thought  that  some  contaminant  might 
exist  which  could  affect  the  diffusion  process.  Experiments  were 
performed  with  commercial  high  purity  cesium  nitrate  from  Kawecki 
Berylco,  with  ultra-pure  cesium  nitrate  from  Merck  and  with  cesium 
nitrate  purified  in  our  own  laboratory.  Figures  9 and  10  show 
almost  identical  profiles  which  were  obtained  from  Merck  and 
Kawecki  Berylco  cesium  nitrate  respectively.  The  dopants  were 
very  different  in  purity,  the  Merck  having  measured  transitior 
metal  impurities  in  the  ppb  range  and  Kawecki  Berylco  in  the  ppm 
range.  It  is  possible  that  all  cesium  samples  tested  had  some 
unmeasured  impurity  such  as  silica  but  the  results  of  many  experi- 
ments did  not  indicate  that  the  source  of  dopant  had  any  effect 
on  the  final  profile. 


F.  Homogeneity 

In  the  early  experiments  the  optical  quality  and  thus  the 
compositional  homogeneity  were  not  very  good.  In  particular,  vi- 
sible striae  were  present  in  the  rods  due  to  incomplete  mixing  and 
due  to  the  folding-in  of  alkali-poor  surface  layers  as  the  batch 
glass  was  stirred.  Since  the  striae  represented  a slightly  dif- 
ferent composition  from  the  average,  the  pores  in  that  region 
would  be  a different  size.  This  would  affect  both  the  diffusion 
and  the  amount  of  dopant  deposited  in  these  pores.  Interferograms 
showed  that  these  striae  persisted  throughout  the  process  and  could 
be  observed  in  the  final  glass.  - 

Stuffing  experiments  were  performed  on  glass  having  both  good 
and  poor  homogeneity.  The  result  was  that  although  the  striae 
affected  the  local  index,  they  had  very  little  influence  on  the 
general  shape  of  the  profile.  Figure  5 has  a good  profile  but 
the  interferogram  in  Figure  11  shows  the  presence  of  striae  which 
distorts  the  interference  pattern.  Figure  12  , however,  shows 

a poor  index  profile  while  the  interferogram  is  relatively  free  of 
striae  as  seen  in  Figure  13  . Homogeneity  of  the  base  glass 

does  not  appear  to  be  the  controlling  factor  in  profile  control. 


The  result  of  the  stuffing  experiments  has  been  that  it  is 
possible  to  obtain  a very  nearly  parabolic  profile  by  molecular 
stuffing  as  predicted.  It  is  clear  however  that  some  parameter 
or  set  of  parameters  not  included  in  the  theory  nay  have  a strong 
influence  on  the  profile.  These  parameters  do  not  appear  to  be 
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among  the  parameters  we  have  discussed  above.  One  of  the  few 
clues  to  come  out  of  this  work  was  the  observation  that  the  pro- 
files produced  using  certain  batches  were  consistently  better  them 
those  produced  from  other  batches.  Batches  having  almost  identical 
analyzed  composition,  coexistence  temperature  and  other  physical 
parameters  simply  did  not  have  the  same  stuffing  characteristics. 

At  the  present  time  we  do  not  have  an  answer  to  the  problems  of 
controlling  the  profile  adequately  and  consistently.  Work  is  being 
continued  in  this  area  with  support  from  private  industry. 

V Profile  Measurements 

A.  Refractive  Index  Profile  Measurements 

Because  of  differing  requirements  at  several  stages  of  this 
project  it  was  necessary  to  develop  three  techniques  for  measuring 
index  profiles  in  fibers  and  in  preforms.  Crude  but  quick  measure- 
ments can  be  made  on  the  preforms-by  measuring  the  angular  deviation 
of  a laser  beam  passing  through  a prism  made  from  the  preform.  More 
accurate  measurements  can  be  made  on  the  preforms  using  a Michelson 
interferometer  operated  in  the  Twyman-Green  mode,  i.e.,  parallel 
illumination.  The  most  accurate  way  of  measuring  the  index  pro- 
files in  fibers  is  the  interference  method  using  an  interference 
microscope.  Because  of  limited  resources,  however,  we  chose  the 
much  simpler  and  less  costly  technique  of  measuring  the  near  field 
intensity  distribution  which  is  related  to  the  index  profile. 

In  our  initial  experiments  we  were  interested  in  obtaining 
quick,  moderately  accurate  measurements  of  the  preform  profiles. 

A simple  way  of  doing  this  is  to  measure  the  deviation  of  a laser 
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beam  passing  through  the  apex  of  a prism  made  from  the  preform. 

The  prism  geometry  is  shown  in  Figure  14  and  a schematic  of  the 
measuring  apparatus  is  shown  in  Figure  15. 

The  procedure  for  sample  preparation  and  measurement  is 
relatively  simple.  A 90°  prism  angle  is  ground  from  the  cylindri- 
cal preform  with  the  angle  bisected  by  the  axis  of  the  cylinder. 

A surface  obtained  by  a final  grind  with  #600  grit  metallurgical 
grinding  paper  is  sufficient.  The  sample  is  placed  with  its  axis 
vertical  in  the  sample  cell  which  has  been  filled  with  paraffin  oil 
(n  = 1.471).  A weakjy  focused  laser  beam  (diameter  =.l  mm)  passes 
partly  through  the  apex  of  the  prism  and  partly  through  the  paraf- 
fin oil  alone.  A zero  reading  is  obtained  on  the  micrometer  by 
rotating  the  mirror  until  the  beam  passing  through  the  oil  alone 
is  reflected  back  on  itself  as  observed  by  the  spot  appearing  on 
the  focusing  lens.  The  beam  passing  through  the  prism  is  then  made 
to  retrace  its  path  by  further  rotating  the  mirror  clockwise  or 
counter-clockwise  depending  on  whether  the  glass  index  is  above 
or  below  that  of  the  paraffin  oil.  This  procedure  is  repeated 
for  various  positions  on  the  radius  of  the  preform  to  obtain  the 
deviation  angle  6 as  a function  of  radial  position. 

The  prism  deviation  angle  can  be  related  to  the  refractive 
index  of  the  prism  and  the  surrounding  medium  vising  geometrical 
optics.  Solving  for  6 one  obtains: 

6 = <J>.  - ot  + sin  ^ {sin [a  - sin”^(— ~ sin  $.)]} 

i n n».  i 

where  the  angles  are  defined  in  Figure  14  and  n and  n'  are  the 

index  of  paraffin  oil  and  the  glass  respectively.  For  the  case 
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of  a 90°  prism  angle  and  a 45°  angle  of  incidence  we,  have  the 
following  simple  relation: 

n'  = n /sin1  (5  + ^-)  + .5 

The  accuracy  of  these  measurements  was  estimated  to  be  better 
than  ± 0.0005  for  an  index  near  1.47,  the  chief  source  of  error 
being  the  accuracy  with  which  one  could  judge  the  moment  when  the 
laser  beam  was  retracing  its  path.  This  estimate  of  the  accuracy 
was  confirmed  by  comparison  of  index  measurements  made  on  our 
samples  in  our  and  in  another  laboratory  (American  Cystoscope 
Makers,  Inc.,  of  Stamford,  Connecticut). 

For  making  precise  observations  of  the  variation  in  index  in 
a cross-section  of  a preform, a Twyman-Green  interferometer  can  be 
used.  A schematic  of  the  interferometer  is  shown  in  Figure  16. 

The  plates  of  the  interferometer  are  adjusted  to  give  an  inter- 
ference pattern  in  the  viewing  plane.  The  sample  is  placed  in  the 
path  between  the  beamsplitter  and  one  of  the  mirrors.  The  observed 
deviations  in  the  interference  pattern  are  due  to  the  additional 
optical  pathlength  in  the  sample.  A sample  having  accurately 
ground  and  polished  flat  faces  and  a uniform  index  will  give  an 
interference  pattern  as  shown  in  Figure  17  which  is  the  photo- 
graph of  the  interference  pattern  from  a homogeneous  piece  of  the 
-untreated  base  glass.  Figures  11  and  13  show  the  interference 
patterns  obtained  from  graded  index  Molecular  Stuffed  samples. 

The  first  of  these  shows  a great  deal  of  inhomogeneity  in  the 
index  profile- due  to  striae  existing  in  the  base  glass  prior  to 
processing.  The  second  of  these  shows  a much  more  uniform  profile 
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but  there  is  still  some  evidence  of  inhcmogeneity  near  the  center 
of  the  sample.  Although  quantitative  results  can  be  obtained  from 
the  interferograms,  we  have  only  used  them  for  qualitative  obser- 
vations to  date. 

B.  Fiber  Index  Profile  Measurements 

Because  of  the  small  diameter  of  the  fibers  (175u  or  less) 
it  is  difficult  to  make  direct  index  measurements.  Although  in- 
terference microscopes  exist  which  will  measure  index  variations 
in  small  samples,  these  instruments  are  costly  and  there  were  not 
any  available  to  our.  laboratory.  An  indirect  technique  of  measure- 
ment is  possible  which  consists  of  measuring  the  near  field  inten- 
sity profile  in  a fiber  for  which  all  modes  are  equally  excited 
with  the  exception  of  leaky  modes.  It  has  been  shown" that  the 
near  field  intensity  profile  is  related  to  the  index  profile 
according  to  the  relation: 

P (r)  = n2 (r)  - n2(2) 

P(o)  n2 (0)  - n2 (2) 

where  P(o)  and  P(r)  is  the  intensity  at  the  fiber  center  and  at 
position  r respectively;  n(0),  n(2)  and  n(r)  are  the  indexes  at 
the  center,  in  the  clad,  and  at  position  r respectively. 

The  apparatus  for  measuring  the  near  field  intensity  pro- 
file of  a fiber  is  shown  schematically  in  Figure  18.  A measure- 
ment  is  made  simply  by  exciting  the  modes  of  a fiber  with  an  in- 
coherent source,  projecting  the  image  of  the  flat  output  end  of 
the  fiber  and  scanning  in  the  image  plane  with  a detector  having 


a small  aperture.  The  difficulty  in  such  a measurement  is  insur- 
ing that  all  ordinary  modes  are  equally  excited  and  that  the  leaky 
modes  are  not  present.  By  making  measurements  using  fibers  of  dif- 
ferent lengths,  one  can  determine  a length  suitable  for  attenuating 
the  lossy,  leaky  modes  while  not  losing  appreciable  power  from  the 
propagating,  higher  order  modes.  For  our  fibers,  one  or  two  meters 
is  usually  sufficient  for  measurment.  Figure  19  shows  index 
profiles  measured  on  a preform  and  a fiber  drawn  from  that  preform. 
The  observed  differences  may  be  real  or  could  be  due  to  errors  in 
the  preform  measurement.  The  steep  portion  of  the  step  profile 
is  difficult  to  measure  accurately  with  the  prism  method.  It  would 
be  useful,  however,  to  confirm  the  results  by  comparison  with 
results  from  an  interference  measurement. 

VI . Physical  Property  Measurement  for  Unstuffing  Bath  Concentration 
Control 

The  computer  studies  of  the  stuffing  procedure  indicated  that 
a few  step  changes  in  stuffing  concentration  with  time  would  be 
sufficient  to  acheive  a parabolic  index  profile  with  a standard 
deviation  of  less  than  .1%  in  index.  The  saturation  concentrations 
of  cesium  nitrate  in  several  solvents  as  a function  of  temperature 
give  a sufficient  range  of  concentrations  to  obtain  the  desired 
profile  control.  These  solubilities  were  measured  for  a number 
of  solvents  and  the  results  are  shown  in  Figure  20.  By  control- 
ling the  temperature  of  the  dopant  solution  the  concentration  is 
conveniently  controlled. 
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vn.  Fiber  Propagation  Characteristics 
A.  Rayleigh  Scattering 
An  upper  limit  to  the  Rayleigh  scattering  in  both  bulk  glass 
and  fibers  may  be  obtained  by  performing  a Landau-Placzek  scat- 
tering measurement  at  a scattering  angle  of  90°.  This  technique, 
combined  with  a measurement  of  total  scattering  and  the  angular 
scattering  distribution,  is  useful  in  separating  contributions  to 
scattering  from  Rayleigh  scatterers,  which  are  related  to  intrin- 
sic bulk  properties,  and  scattering  from  core  diameter  perturba- 
tions and/or  needle-like  scattering  particles.  Since  most  of  the 
Rayleigh  scattering  is  related  to  intrinsic  bulk  properties  of  the 
glass,  it  represents  a lower  limit  to  the  attenuation. 

The  measurement  of  Landau-Placzek  ratios  in  glasses  has  been 
described  thoroughly  in  the  literature.  ' Briefly,  the 

light  scattered  by  the  glass  from  a collimated  laser  beam  is  col- 
lected at  a convenient  scattering  angle  such  as  90°.  The  light 
is  frequency  analyzed  by  a scanning  Fabry-Perot  interferometer 
which  resolves  the  unshifted  Rayleigh  line  and  the  Brillouin 
lines  which  are  shifted  upward  and  downward  in  frequency  due  to 
scattering  from  acoustic  waves.  A schematic  of  a typical  experi- 
mental setup  is  shown  in  Figure  21.  . The  Brillouin  intensities 
are  a weak  function  of  the  glass  composition,  while  the  Rayleigh 
scattering  is  a function  of  both  composition  and  thermal  history 
of  the  glass.  The  absolute  Brillouin  intensity  can  be  calculated 
if  the  stress  optical  coefficients  are  known  for  the  glass  of 
interest.  The  Brillouin  intensities  can  then  be  used  as  a cali- 
bration for  estimating  the  Rayleigh  intensity.  If  the  coefficients 
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are  not  known,  a less  accurate  result  may  be  obtained  by  using 
a measurement  of  a known  glass  such  as  fused  silica  to  calibrate 
the  sample  measurement.  The  latter  technique  was  used  to  measure 
the  upper  limit  of  the  intrinsic  scattering  in  both  preform  and 
fiber  samples  of  Molecularly  Stuffed  glass. 

The  scattering  attenuation  is  given  by 


(1) (2)  •lB 
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(1)  (2) 

where  a , a , are  the  scattering  attenuation  of  the  sample 

9 0 " • 

and  of  fused  silica  respectively;  IB(1*,  IB^2)  are  the  Brillouin 
intensities  of  the  sample  and  fused  silica  and  are 


the  Landau-Placzek  ratios  of  the  sample  and  fused  silica  respec- 
tively where  the  Landau-Placzek  ratio  is  defined  by  = IR/2lfi. 

(2) 

a is  known  and  is  approximately  1 dB/km  at  0.85y. 

Measurements  performed  on  Suprasil  II  fused  silica  and 
cesium  doped  Molecular  Stuffed  preforms  indicated  a scattering 
loss  for  the  Molecular  Stuffed  preform  of  approximately  2 dB/km 
at  .85y.  Measurements  performed  on  fibers  involve  a slight  com- 
plication due  to  the  fact  that  the  large  numerical  aperture  of  the 
fibers  means  that  light  collected  at  90®  from  the  fiber  axis  and 
beam  direction  may  actually  be  scattered  at  a range  of  angles. 

This  is  due  to  the  fact  that  if  all  of  the  fiber  nodes  are  ex- 


cited, then  light  is  propagating  in  the  fiber  at  all  angles  al- 
lowed within  the  NA  of  the  fibers.  This  results  in  a broadening 


of  the  Brillouin  peaks  and  one  must  be  careful  to  use  the  inte- 
grated Brillouin  intensities  in  calculating  the  Landau-Placzek 
ratios.  Also, because  of  the  small  size  of  the  fiber,  it  is  dif- 
ficult to  avoid  collecting  light  scattered  from  the  fiber  surface 
This  may  be  minimized  by  carefully  stripping  the  cladding  modes 
from  the  fiber  with  a suitable  mode  stripper.  Measurements  made 
on  a Molecularly  Stuffed  fiber  indicated  a loss  of  approximately 
2.5  dB/km  which  is  in  reasonable  agreement  with  the  bulk  measure- 
ment considering  the  uncertainties  in  the  measurement. 

B.  Total  Light  Scattering  Measurements 

In  addition  to  Rayleigh  scattering,  fibers  may  exhibit 
significant  scattering  due  to  a variety  of  particles  or  local 
dielectric  variations  present  in  the  fiber.  In  this  case  a scat- 
tering measurement  at  one  scattering  angle  is  not  sufficient 
to  obtain  the  total  scattering  at  all  angles.  An  integrating 

sphere  or  similar  device  is  required  to  make  such  a measurement. 

. . 4 

A simple  integrating  cube  was  described  by  Tynes  that  is 

very  convenient  for  measuring  the  total  scattering  in  fibers. 

This  integrating  cube  is  constructed  of  six  matched  silicon 
solar  cells.  Two  opposing  faces  of  the  cube  have  small  holes 
to  allow  the  entrance  and  exit  of  the  fiber.  The  cell  is  filled 
with  an  index  matching  fluid  such  as  glycerin  that  couples  out 
the  light  scattered  from  the  core  into  the  clad  in  the  near  for- 
ward direction.  It  is  thus  possible  to  collect  nearly  all  of  the 
light  scattered  from  the  core  with  the  exception  of  that  striking 
small  dead  areas  in  the  cube  near  the  edges  and  near  the  fiber 
entrance  and  exit  holes. 
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A schematic  of  the  experimental  setup  is  shown  in  Figure  22. 
The  light  from  a laser  or  other  bright  light  source  is  focused 
on  the  end  of  the  fiber.  Before  entering  the  scattering  cell, 
the  fiber  is  passed  through  a tube  filled  with  glycerin  to  strip 
out  accumulated  cladding  modes.  The  fiber  is  passed  through  the 
cell  which  has  a length  of  1 cm.  The  light  scattered  b^  the  core 
in  that  one  cm  is  detected  and  the  resultant  signal  is  read  from 
a voltmeter.  The  scattering  as  a function  of  position  in  the 
fiber  is  measured.  Regions  of  minimum  scattering  for  a good  fiber 
are  expected  to  be  representative  of  the  intrinsic  Rayleigh  scat- 
tering while  high  scattering  regions  have  large  contributions  from 
other  scattering  sources.  The  absolute  scattering  attenuation  at 
any  point  may  be  determined  by  first  measuring  the  scattering  in- 
tensity and  then  breaking  the  fiber  and  using  the  cell  to  measure 
the  light  being  transmitted  in  the  fiber  at  that  point.  The  scat- 
tering attenuation  for  small  scattering  is  then  given  simply  by 

C I«. 

ag (dB/km)  = 10°  1°51oe 

o 

where  I and  I are  the  scattered  and  transmitted  intensities 
s o 

respectively. 

As  a check  on  the  absolute  magnitude  of  the  measured  scatter- 
ing intensities  a silicone  clad  Suprasil  II  fiber  with  total  at- 
tenuation at  0.85u  of  approximately  20  dB/km  was  measured.  The 
average  scattering  loss  was  approximately  6 dB/km  which  is  the 
order  of  magnitude  expected  for  such  a fiber.  The  minimum  scat- 
tering loss  for  this  fiber  was  about  3 dB/km  which  is  higher  than 
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the  expected  intrinsic  Rayleigh  scattering  loss  but  probably  con- 
tains excess  contributions  from  the  silicone  silica  interface 
because  of  the  large  NA  used  in  making  this  scattering  measurement. 

A mdlecular  stuffed  fiber  with  total  attenuation  at  0.85u  of 
approximately  25  dB/km  when  measured  with  light  launched  at  0.07  NA 
was  examined  to  determine  the  minimum  scattering.  Five  one-cm-long 
regions  having  an  average  scattering  attenuation  of  2.7  dB/km 
measured  with  a launching  NA  of  0.09  were  found  by  examining  five 
short  sections  of  fiber  chosen  at  random  from  the  20-meter-long 
fiber.  A second  measurement  was  made  at  the  same  points  using  a 
launching  NA  of  0.1R.  These  measurements  indicated  an  increase 
in  average  scattering  to  a value  of  4.2  dB/km.  Tiis  is  probably 
due  to  increased  scattering  from  the  core  clad  interface.  These 
low  scattering  regions  were  further  analyzed  with  respect  to  the 
angular  distribution  of  the  scattered  light  as  is  described  in  the 
next  section.  The  angular  distribution  found  was  identical  to 
that  expected  for  Rayleigh  scattering.  The  magnitude  of  the  to- 
tal scattering  is  consistent  with  that  found  in  tie  Rayleigh- 
Brillouin  measurements  reinforcing  the  conclusion  that  the  intrin- 
sic scattering  for  these  fibers  is  less  than  2.5  dB/km  at  0.85u. 

C.  Angular  Distribution  of  Scattered  Light 

The  measurement  of  the  angular  distribution  of  the  scattered 

light  can  be  used  to  obtain  information  about  the  nature  of  the 

scattering  mechanism.  Rayleigh  scattering  from  a linearly  pola- 

2 

rized  beam  is  characterized  by  a (1  + cos  0)  angular  distribution 
factor.  Both  core  diameter  variations  and  long  dielectric  needles 
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can  give  rise  to  a strong  forward  scattering  contribution.  Measure- 
ments of  the  angular  scattering  for  fibers  was  reported  by  Rawson  ^ 
along  with  his  interpretation  identifying  both  the  Rayleigh  con- 
tribution and  contributions  from  dielectric  cylinders. 

The  angular  distribution  of  scattered  light  from  molecular 
stuffed  fibers  was  measured.  Particular  interest  was  paid  to 
regions  identified  in  the  previous  total  scattering  measurements 
as  strong  scattering  regions  and  weak  scattering  regions.  A 
schematic  of  the  experimental  setup  is  shown  in  Figure  23.  The 
light  from  a laser  was  focused  on  the  end  of  a fiber.  The  fiber 
was  passed  through  a_ cell  filled  with  paraffin  oil.  The  cell  was 
centered  on  a rotatable  table  whose  angular  orientation  could  be 
read  to  one  minute  of  arc.  Black  paint  applied  tc  all  but  the 
central  one  cm  of  the  fiber  served  both  to  mask  all  but  that 
section  from  the  detector  optics  and  stripped  accumulated  cladding 
modes.  By  measuring  scattered  intensity  as  a function  of  the  orien- 
tation of  the  fiber  with  respect  to  the  axis  of  the  collection 
optics,  and  making  the  appropriate  corrections  for  changes  in  the 
collection  volume  with  orientation,  the  angular  intensity  distri- 
bution could  be  obtained. 

2 

It  was  found  that  the  weak  scattering  regions  had  a (1+cos  Q) 
angular  distribution  characteristic  of  Rayleigh  scattering.  A 
typical  intensity-versus -angle  plot  for  one  of  those  regions  is 
shown  in  Figure  24.  The  strong  scattering  regions  however  had 
both  the  Rayleigh  scattering  contributions  as  well  as  a strong  for- 
ward scattering  contribution  which  could  be  orders  of  magnitude 
greater  than  the  Rayleigh  contribution.  A typical  angular  dis- 


V 


24. 

tribution  for  one  of  the  strong  scattering  regions  is  shown  in 
Figures  25  and  26. 

If  the  source  of  the  strong  forward  scattering  is  assumed  to 
be  dielectric  cylinders  as  in  the  interpretation  of  Rawson,  5 one 
can  estimate  the  length  of  the  cylinders  from  the  angular  width  of 
the  forward  scattering  peak  according  to  the  expression 
el/2  = t-8856  x /anQn1/2 

where  0^is  the  angular  width  at  1/2  maximum,  X is. 5145  y , nQ  is 
1.47.  The  index  of  paraffin  oil  and  L is  the  length  at  the  scat- 
tering cylinders.  For  the  scattering  region  of  Figure  25  the 
length  of  the  cylinders  is  estimated  to  be  approximately  250y.  If 
the  cylinders  are  assumed  to  be  the  result  of  the  stretching  dur- 
ing fiber  drawing  of  small  spherical  inhomogeneities  existing  in 
the  preforms,  it  is  estimated  that  the  diameter  of  the  spheres 
in  the  preforms  would  be  approximately  .25y.  An  examination  of 
the  preforms  in  an  electron'  microscope  revealed  the  presence'  of 
spheres  approximately  .25y  in  diameter.  The  source  of  these 
spheres  is  not  known  but  may  be  due  to  silica  gel  formed  during 
the  molecular  stuffing  process  and  not  fully  dissolved  in  the 

sintering  step  or  in  fiber  drawing. 

D.  Pulse  Dispersion  Measurements 

Due  to  the  fact  that  the  various  modes  in  a multimode  fiber 
do  not  travel  down  the  fiber  with  equal  velocity,  the  phenomenon 
of  pulse  broadening  is  observed.  If  one  considers  the  modes  to 
be  represented  by  rays  launched  into  the  fiber  at  specific  angles 
to  the  fiber  axis,  it  is  easy  to  see  the  source  of  the  pulse 
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broadening.  Simply  stated,  a ray  traveling  at  a small  angle  with 
respect  to  the  fiber  axis  travels  a shorter  distance  than  one 
traveling  at  a steep  angle.  The  resultant  pulse  broadening  or 
dispersion  can  be  measured  directly  by  launching  a narrow  pulse 
into  a fiber  and  measuring  the  width  of  the  pulse  after  passing 
through  a length  of  fiber. 

Apparatus  for  measuring  pulse  dispersion  was  built  in  our 
laboratory  with  a design  similar  to  that  described  by  Andrews.6 
A schematic  of  the  measurement  apparatus  is  shown  in  Figure  27. 

The  light  source  is  an  injection-  GaAs  laser  (RCA  SG-2001)  driven 
by  5 ampere  5 nsecond  pulses  at  a repetition  rate  of  2 to  5 kHz. 

The  light  pulse  (A=.9y)  is  collected  by  a 20X  microscope  objec- 
tive and  coupled  into  the  fiber  with  a 100X  objective.  The  light 
exiting  the  fiber  is  detected  by  either  a PIN  or  an  avalanche 
silicon  photodetector  whose  output  is  recorded  by  a sampling 
oscilloscope.  The  triggering  of  the  oscilloscope  is  delayed  with 
respect  to  the  launched  pulse  to  allow  for  the  transit  time  of  the 
pulse  through  the  fiber. 

Dispersions  measurements  were  made  on  a 1 km-long  step  index 
Molecular  Stuffed  fiber  having  an  attenuation  at  . 9y  of  approxi- 
-mately  20  dB/km.  Figure  28  shows  the  pulse  shape  after  traver- 
ing  several  meters  and  after  traversing  1 km  of  fiber.  The  pulse 
from  the  short  fiber  has  a width  at  half  maximum  of  .8  nsec  which 
represents  the  instrumental  width  of  the  system.  The  pulse  from 
the  long  fiber  is  highly  asymmetric  due  the  the  high  attenuation 
of  the  higher  order  modes  and  has  a width  of  13  nsec.  This  im- 
plies a bandwidth  for  the  fiber  of  the  order  of  70  to  80  Mbits/sec-km. 


Molecular  Stuffed  graded  index  preforms  of  sufficient 
quality  and  size  have  not  been  available  for  drawing  fibers 
suitable  for  dispersion  measurements.  It  would  be  expected  that 
these  preforms  would  give  improved  bandpass  over  the  step  pro- 
filed preforms. 
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CONCLUSION 

The  purpose  of  this  research  was  the  development  of  the 
Molecular  Stuffing  process  to  produce  low  loss  optical  fibers 
having  a parabolic  index  profile.  The  first  stage  of  the  research 
was  the  development  of  a diffusion  theory  applicable  to  a porous 
cylinder  immersed  in  a dopant  solution.  The  second  stage  involved 
the  determination  of  diffusion  coefficients  solubilities  and 
related  physical  parameters  needed  for  application  of  the  theory. 
The  final  stage  was  the  production  of  preforms  having  parabolic 
index  profiles. 

The  development  of  the  diffusion  theory  was  very  successful. 

It  was  discovered  that  by  a series  of  step  changes  in  dopant 
concentration  with  time  or  by  a step  change  followed  by  a 
ramp  in  concentration  one  could  theoretically  obtain  a parabolic 
or  if  desired,  a near-parabolic  profile  to  within  a standard 
deviation  of  10”®.  This  was  verified  by  direct  computer  cal- 
culation. 

As  a dopant,  CSNO3  in  aqueous  solution  proved  to  have  very 
desirable  characteristics.  It  could  be  purified  relatively  easily. 

It  was  found  that  solution  concentration  could  be  well  controlled 
by  controlling  the  temperature  of  a saturated  solution. 
Precipitation  and  controlled  disssolution  of  the  CSNO3  could  be 
obtained  by  temperature  change  and  the  use  of  organic  solvents. 
CSNO3,  which  becomes  CS2O  in  the  final  glass  also  has  the  index 
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modifying  effects  necessary  for  use  in  optical  fibers.  In  short, 
it  met  all  the  requirements  for  use  in  Molecular  Stuffing. 

A large  number  of  preforms  were  produced  with  varying 
success.  The  fact  that  several  profiles  were  produced  with  nearly 
parabolic  profiles  has  convinced  us  that  Molecular  Stuffing  as 
it  has  been  developed  can  be  successful.  There  are  however, 
important  effects  which  are  not  accounted  for  in  our  theories 
and  further  work  is  required  to  control  these  effects.  This 
work  is  continuing  with  support  from  private  industry  and  it  is 
expected  that  the  results  will  be  positive. 
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Pig.  1.  Plow  diagram  of  the  Molecular  Stuffing  process. 
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Fig.  2 An  early  graded  index  profile  showing  various  profile 
errors  including  a flat  top,  high  index  at  the  edges 
and  a skewed  profile  shape.  Rod  2097. 
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Pig.  5 The  profile  of  a rod  given  an  HF  wash.  The  standard 
derivation  of  the  data  from  the  parabola  shown  by  the 
solid  line  is  4%  for  the  central  75%  of  the  profile. 
Rod  2518. 
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Fig.  9, 


The  profile  of  a rod  doped  with  Kawecdti-Berylco  CSNO3. 
The  profile  is  moderately  good  in  shape  hut  the  edges 
are  too  high.  Rod  2892. 
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Fig.  10.  The  profile  of  a rod  doped  with  ultra-pure  CSNO3  from 
Merck.  The  profile  is  nearly  identical  to  that  shown 
in  Fig.  9.  Rod  2889. 
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Fig.  11  Interferogram  of  a rod  cross-sectiom  havirg  a parabolic 
profile.  General  inhomogeneity  in  the  base  glass  is 
observed. 
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Fig.  12.  The  profile  of  a preform  prepared  from  a homogeneous 
glass  rod.  The  profile  has  a flat  top  and  high  edges 


Fig.  14.  PRISM  GEOMETRY  FOR  INDEX  MEASUREMENT 


Figure  15:  Index  Profile  Measuring  Apparatus 

A.  Sample  Prism 

B.  Mirror 

C.  Optical  Cell 

D.  Matching  Fluid  - Paraffin  Oil 

E.  Ball  bearing,  « 1/2  diameter 

F.  Micrometer 


Figure  16.  Michelson  Interferometer 
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Figure  19. 

The  solid  line  is  the  plot  of  the  index  profile  of  a fiber 
measured  using  the  near  field  technique.  The  triangular  data 
points  show  the  profile  of  the  rod  from  which  the  fiber  was 
drawn. 
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Figure  22, 


Schematic  of  the  total  scattering  measuring 
apparatus . 


Figure  24.  The  angular  distribution  of  light  scattered  from  a 
fiber  section  exhibiting  only  the  1 + cos^S  depen- 
pendence  due  to  Rayleigh  Scattering.  The  1 + cos^9 
shape  is  distorted  due  to  refraction  at  the  core 
clad  and  clad  oil  interfaces.  The  noisy  signal  in 
the  near  forward  direction  is  due  to  the  difficulty 
of  collecting  the  light  scattered  in  the  forward 


direction. 
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Figure  25.  The  forward  scattering  peak  of  a strong  scattering 


region  of  a fiber.  The  intensity  scattered  into 
the  forward  peak  is  21  times  that  due  to  the  inte 
grated  Rayleigh  scattering  alone. 
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Figure  26.  The  angular  distribution  of  light  scattered  from  a 
strong  scattering  region  in  a fiber.  The  forward 
peak  is  offscale  in  this  plot. 
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Figure  27.  Schematic  of  the  Dispersion  Measuring  Apparatus 
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